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ON TRANSONIC FLOW PAST A WAVE-SHAPED WALL ’ 

By Cabi, Kaplan 


SUMMARY 

The present report is an extension of a previous investigation 
(described in NA.CA Rep. 1069) concerned with the solution 
of the nonlinear differential equation for transonic flow past 
a wavy wall. In the present work several new notions are 
introduced which permit the solution of the recursion formulas ■. 
arising from the method of integration in series. In addition, 
a novel numerical test of convergence, applied to the power ' 
series (in the transonic similarity parameter) representing the 
local Mach number distribution at the boundary, indicates 
that smooth symmetrical potential flow past the wavy wall is ' 
no longer possible once the critical value of the stream Mach 
number has been exceeded. 


INTRODUCTION 

In NACA Report 1069 (ref. 1) the writer considered the 
problem of two-dimensional transonic flow past an infinitely 
long sinusoidal solid boundary. The problem was treated 
in the physical plane and the purpose was to investigate 
whether or not the flow past the wavy wall remains a smooth 
symmetrical type of potential flow when the undisturbed- 
stream Mach number exceeds its critical value. By a 
smooth type of potential flow is meant one for which the 
velocity potential, say, and its first derivatives are single- 
valued and continuous; that is, there are no discontinuities 
of the nature of shock waves. 

Initially, the Brand tl-Busemann small-perturbation method 
was applied and the velocity potential developed inclusive 

of the third power in the “thickness” coefficient « =- ^ 

A/^TT 


where 

a amplitude of wavy wall 

X wave length of wavy wall 

The velocity potential was then referred to the critical speed 


of sound CcT, the coefficient e was replaced by 


k(l-My^ 

T+ 1 


where k = 


(7 + l)e 
(l-My^ 


is the transonic similarity parameter, and 


terms involving powers of 1 —MJ higher than the first were 
neglected (main assumption for transonic flow) . This simpli- 
fied or transonic form of the Brand tl-Busemann solution was 
shown to be identical with the one obtained by solving the sim- 
plified nonlinear differential equation for transonic flow, with 
the boundary condition taken not at the wave-shaped wall but 
at the flat plate corresponding to vanishing amplitude. The 
calculation was carried through the sixth power in the tran- 


sonic similarity parameter k and corresponds to the in- 
superable task of obtaining the Brandtl-Busemann solution 
to the sixth power in the thickness coefficient e. Thus each 
iteration step of the Brandtl-Busemann method contributes 
to the transonic form of the solution, which may therefore be 
considered a result of thin-profile theory with disturbances 
not necessarily small compared with l—Mff. The main 
conclusion reached in reference 1 was that the transonic 

4 

similarity parameter k must be less than ^ — a value still 

o 

somewhat greater than the critical value 0.83770 calculated 
there. 

The purpose of the present work is to express the solution 
of the problem of transonic flow past the wavy wall in a 
form more suitable for general considerations and to prove 
that the assumed smooth symmetrical type of potential flow 
cannot exist at stream Mach numbers beyond the critical 
value. 


SYMBOLS 


^,y 

a 

X 

X/27T 


7 

u 

M 


k 

^cr 

<t> 

/(*, y) 
ff{y) 

A" I 


nondimensional I’ectangular Cartesian coordinates 
amplitude of wavy wall 
wave length of wavy wall 
reference element of length 


‘thickness” coefficient of wavy wall. 


X/27t 


ratio of specific heats at constant pressure and con- 
stant volume 
stream velocity 
local Mach number 
undisturbed-stream Mach number 

(7+l)« 


transonic similarity parameter. 


(1-MT' 


critical speed of sound 
velocity potential of flow 
transonic disturbance potential 
functions of y only, related to f(x,y) 
numerical coefficients 


CP 

generating functions of fc, u+r 

r=0 


Atop functions of dummy variable r, 23 lower 

n = l,7n 

label starts from 1 when m is negative and from 
m when m is positive 

Brimes denote differentiation with respect to independent 
variable. 


* Supersedes NACA TN 2748, ‘'On Transonic Flow Past a Wave-Shaped Wall” by Carl Kaplan, 1952. 


279836—54 
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ANALYSIS 

GENERAI. FORMULAS 

If the undisturbed stream is in the direction of the positive 
a:-axis, then the velocity potential <l> referred to the critical 
velocity Ccr can be written as (see ref. 1) 

where the second term on the right-hand side is a disturbance 
velocity potential and implies that terms involving powers of 
1 ~MJ higher than the first have been neglected. The differ- 
ential equation for transonic flow satisfled by the function 
J{x,y) is obtaiped from the general differential equation for 
compressible flow and takes the following simplified non- 
linear form: 

dy' da: da;' ^ ’ 


the following system of second-order nonlinear ordinary 
differential equations for/„ results: 

J n-l 

jn"—n^Jn= — y wS w (n— — 

4 m=l 

(w=l,2, . . . oo) (4) 

" m=l 

Before proceeding to the solution of these equations, several 
formulas of general interest and subsequent use are given. 
They have been derived in reference 1. The local Mach 
number distribution is given by 

(5) 

The equation from which the critical value of the transonic 
similarity parameter is calculated follows from equation (5) 
by taking a;=0, j/=0, and M=l; that is, 


The boundary conditions to be fulfilled by f(x,y) are as 
follows : 


I- i 

_ h 


(at y= oo) 


oo <(a;<( oo) 


( 2 ) 


Here x and y are nondimensional rectangular Cartesian 
coordinates simply related to the physical plane coordinates 
X and Y by means of the transformation 

x^X 



y-O 


or with the aid of equation (3), 

Sn/„(0) = 1 

71<=1 


The pressure coefficient 


is given by 


Op. 


Ma> 


P—P<. 

\p.U‘ 


Op 


Ma, 


2 

T+1 




( 6 ) 


( 7 ) 


y = (l-M„T'F 

and the equation of the wavy wall is Y=e cos X or 
y—{l—Mjy^t cos X. Clearly, j{x,y) involves only the 
variables x and y and the transonic similarity parameter k. 

The most general form for the function j{x,y) to ensure 
symmetrical potential flow past the wavy wall is the following 
(see fig. 1) : 

— sin nx (3) 

71 = 1 

where the are functions of y only. When this form for 
j{x,y) is substituted into equation (1) and the coefficients of 
the separate harmonic terms sin nx are equated to zero. 


Y 



INTEGRATION IN SERIES FOR THE FUNCTIONS /•„ 

In reference 1, equations (4) were solved by an iteration 
procedure and approximate expressions for /i to /e were ob- 
tained. An examination of these expressions showed that 
the general form of /„ is 

71 — 1 

(» 2p+7i-2 CO 

/n=S S (rt=l,2, . . .a,) (8) 

p=0 g=0 T=p 

where, if p = 0, the upper limit of q is n—1 and, if the 

upper limit of q is 2jo+n — 2. The four-labeled coefficients 
Al ?are real numbers calculated from recursion formulas 
obtained from the system of differential equations (4) and 
the boundary condition at the surface of the wavy wall. 
The boundary condition at y= is automatically satisfied 
by the form of ; whereas the boundary condition at the 
wall takes the form 

(:;;)} (9) 

Inserting the expression for f„ given by equation (8) into 
equations (9) yields immediately the following results: 

Aro=i \ 

i:(2p-fn)ylS ?= S A? ? (ti= 1, 2, . . . -) f 

p=0 1 
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where, if n=l, the lower limit of p on the right-hand side 
is unity and, thelowerlimitofpiszero. Also,ifn=l, 

the upper limit r of p goes from 1 to » and, if r goes 

from 0 to 05 . 

By elementary manipulation of series, the second of equa- 
tions (10) can be replaced by the following more useful forms: 


nAz ! 


(n=2, 3, 


nAz°= S ZK22>+«-)^or” o’ ■ 


( 11 ) 


where in the first term on the right-hand side the lower 
label p starts from 1 when n=l and from 0 when n9^\. 

In reference 1, recursion formulas were derived for the 
coefficients . In the present report a much more 
significant approach is introduced. Note that equation (8) 
can be rewritten in the following form : 


Let 

Ai o=S 

n = l 

then 

^io' = S 

71 = 1 

and 

{Ai o')^=S S rn(p,—m)AZ-% AlziA 

71=2 7n=l 


By observation, it can be seen that the right-hand side of 
equation (15) is equal to {rAi o'Y. Equation (15) can thus be 
replaced by the following first-order nonlinear differential 
equation : 

{rA,,'f-%rA,o’ + SA,o=0 (16) 


The solution of this equation is 


71 — 1 

« 2p +77— 2 

fn = 'E>{ke-''y^+” Z) A"/ (12) 

p=fl 9=0 

where 

An V V h2r 

a — A v+r ^ 

r=0 

In a manner similar to that described in reference 1 for the 
coefficients A^ f , recursion formulas can be obtained for the 
power series A"'/. Indeed, the two types of recursion formu- 
las are intimately connected and, for a given value of p, 
the one can be obtained from the other by mere inspection. 
A single recursion formula can be written for the general 
quantity A"/ but the resulting expression is cumbersome 
and serves no practical purpose. It is much more desirable 
to obtain separate recursion formulas for each value of p. 
As examples of procedure, the recursion formulas for p=0, 
1 , and 2 are considered in the sections which follow. 

RECURSION FORMULA FOR 

With p = 0, the recursion formula is (compare with cq. (57) 
of ref. 1) 

2niq+l)AU1={q+l) (9 + 2)6" J + 

7 ft in S" (m-f-l)(ft— m — 1)A"-™-* 

4 9i=0 m=q-~q\ 

(ft=2, 3, . . . oo; ^ = 0, 1, . . . ft— 2) (13) 

where 

3n-2_0 (g = ft-2) 

“ 1 (?7^ft-2) 

This recursion formula can be solved, the solution starting 
with 2=71—2 and descending towards 2=0. Thus, for 
2 =ft— 2, equation (13) becomes 


A^ 0 — 


where c„ is the arbitrary constant of integration. From the 
definition of Ai o it follows that, with r=0. 


and hence 





A\ p' 


(17) 


Note that by definition +"o“ involves the coefficients Al “ 
which are ultimately calculated by means of the auxiliary 
relations, equations (11), engendered by the boundary con- 
ditions at the surface of the wavy wall. In solutions of the 
recursion formulas, therefore, the coefficients +” 0 ° appear as 
undetermined quantities. 


A ' 

The expansion of --+t powers of rA'o" is a nontrivial 
A. 0 

problem which fortunately can be solved with the aid of 
Lagrange’s investigations on the reversion of power series. 
In reference 2, Lagrange’s problem is illustrated by the 
following example: 

To expand 6®“^ in powers of y=xe’’^. The result is expressed 


as 


l+ay-\ 


a(a- 


2! 


2b) , 

y - 


a(g — 3 by 
3! 


y^+- ■ ■+ 


a(a—nb)’^~^ 

ft! 


2/" + - • . 


(18) 


In particular, with a=l, b= — 1, and the solution of 

the equation log J=7/| or J=e‘'« (generally referred to as 
Eisenstein’s problem) is given by the series 


8(ft— 1 ) 24 ;; _?=X1 ("I + 1) (ft— m— l)A’"+j; 

771=0 


(ft=2,3,. . .co) (14) 


Now, multiply both sides of this equation by r” where r is a 
dummy variable and sum from n=2 to n= 00 . Then 

8X)(ft— 1)^" («^ + l) [n—m — l)A"'+^ °A\-jr~} 2 ° 

71 = 1 71 = 2 771 = 0 

( 15 ) 




ft" 


^i(ft— 1)! 




A ' 1 

Then, if f is replaced by .V°o ~ V solution 

of equation (17) is 


re -2 


0 ^ ( A\ 0Nre-:«re-I 
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or 

Hence, by equating coefficients of equal powers of r on both 
sides of this equation, 

(71=2,3, ... co) (19) 

which is the solution of the recursion formula, equation (14). 
The coefficient can be considered as the generating 

function for the set of coefficients A„”^i ?. Thus, equation 
(19) can be written as 




( 20 ) 


Then, by equating coefficients of equal powers of k on both 
sides of this equation and putting HJg=l (see boundary 
condition, eqs. (10)), the following equations are obtained : 


A -I 


0 ^ 

“ n!4”-' 

n -2 


AA-ri= 


71" 


(n— 1)!4 

n' 


J 1 0 

n -1 -^0 1 


(ti— 1)!4 


/d 1 0 

n-l -^0 2 


2!(ti-2)!4’ 


{Al If 


^ ( 21 ) 


Note that from the first of these formulas Al S==- and that 

O 

'from the boundary condition, equations (ll),H^o=j^- 
Consider now q=n — 5. Then equation (13) becomes 
2n(7i-2)H;_“=(7i-l) (ti-2)H:_;+ 


^n'^m(n-m)AZ-lA’^ZZ-1 (22) 

Zi m=2 

Multiply both sides of this equation by r” and sum from 
71=3 to 77= CO. Then 


2 Xj 7 i( 77 — 2)H"_2 7 ’'' = ]2 (ti— 1 ) (71— 2)H;_;r" + 

n= 2 n=l " 


Let 


and 


nr" m (n—m) AZ-l A” zZ-i (23) 

72 = 3 77» = 2 


A,o=i:A:Ar" 

n = l 


A,o=J2AzA,r" 

n = 2 


then 


Hi o'Ma o' = Xl7'”S"»(7*'— 7a + 2)H”_;H" "im" 

71 = 1 7n = l 


It can be shown easily that the last term on the right-hand 

side of equation (23) is equal to ^ r (r-^Hi o'H 2 oO'- Hence, 

equation (23) can be rewritten as the following differential 
equation ; 

2 HH 2 0 "— 2 rH 2 o' = 7-^Hi 0 "— 2?’Hi o' + 2Hi a~\-\r{r^Ai o'H 2 oO' 

( 24 ) 

Now, from equation (16), it follows that 
(rA 

(7-Hi o') 

rA 

Then equation (24) can be written as 


1 0 ^)^ — 8 (t’Hi o' — Hi 0 ) \ 
4 r\/ 4 H 1 0 ^ \ 

^'“^=4=7HT7 { 

io'(rHio')' = 4rH. ) 


(25) 


(^H 2 °? V YrA, o')= 


r(4 — 7 -Hi o') 


V(4-rHi o') 


(7-Hi 0 ')^ 


LetrH2o'=7.4_^^^^, 
iable. Then equation (26) becomes 


6rH, o' + 4H, 0 ] (26) 

7 J where v is the new dependent var- 


V 


A{rAi o') 
or with the aid of equations (25), 


[2r(7-Hi 0 ')'— 6t’Hi o' + 4Hi 0 ] 


Hence 


or 


v'=^ {rA, 0 ')' 


«=| rA, o'+i c, 




4-rA.' 4-M,,' 

The arbitrary constant Ci is determined by differentiating 
this last equation twice and evaluating for r=0. Thus, 


Finally, 


A 2 0 

(H'oT 

H 2 o'=~- [(Hi o')1' + 8 rHi 0 " (27) 


Now, from equation (18) with a— 2, 6= — 1, and ?=c*, it 
follows that 

1 n —3 


?^=2S 


97" 


( 77 — 2 )! 


yn 
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Or with 


o' 

^ AV 


y=4 '’^0° 


and 


Thus, equation (27) becomes 

0 ) n=2 

Hence, replacing by its value from equation (19) and 

equating the coefficients of equal powers of r on both sides 
of this equation yields 


A^A=8 


7)^ -3 1 

ir^(^‘o'’)"-^^V+ 


(m-2)!4 


(AVr 


2 (n-3)!4" 

(71=3,4, ... co) (28) 


the solution of the recursion formula, equation (22). Again, 
equation (28) can be written as 

§»••-. ; ±a; ; r- (±.4! ; k-J -+ 

Then, equating the coefficients of equal powers of k on both 

sides of this equation and putting AJo=l and g==— 

1 o 

yields : 

A n 0 

-2 0 


2(71- 2)!4"-'^2(a- 3)14''-^ 


8a" 


A n 0 A20 J A\ 0 

(7i,-2)!4»-^ " ‘ + (a-3)!4"-‘ ' 

C,jg n — 3 

A n 0 J 2 0 I 

(n-2)!4"-‘ 

qnJl ~Z r 


(30) 


Note that from the first of these formulas A\ S=kv and that 
from the boundary conditions, equations (11), AH—-^- 

I £t 


The procedure followed in order to obtain equation (28) 
for .4^ _S is a general one and with very little difficulty other 
members of the family A",” can be obtained. For example, 


a" 


An 0 4 0 fAlO\n-3A30\ 

-^ln-3 — 0 ) 0 m 

1 


<jn-4 


i ^ji^Mvy-*[iQAv+(AV)r 

From this equation it follows that 


(31) 


A"; 


1 


'3 (a— 3)!4"-‘ ' 2 (a— 4)!4’‘-‘ 


A„U ?=48 


a" 




(a-3)!4"-‘ (a-4)!4 

0 R , 1 

^ 0 1 1_3 (^_4)!^2 (a-5)!J 


?-l- 


M32) 


From the first of these formulas and the boundary conditions, 
7 7 

equations (11), = and A% °= ^g 3 Q- 

At this point, it is noted that the coefficients of the form 
A^ S are calculated from the first formula of each set, equa- 
tions (21), (30), (32), and so forth. A number of this type of 
coefficient have been evaluated (see ref. 1) by means of the 
recursion formulas for the coefficients A" ? themselves. 
They are listed as follows: 


AH — 

1 

0 _ 

7X13 

0 

8 

-^1 0 

72X256 

09 

II 

1 

47 0 

13 

24 

■^X 0 

18X256 

A 4 0 

jH-1 0 — 

7 

48 0 

13X19 

384 

-^1 0 

576X256 

A5 0 

■^1 0 — 

7 



768 




A careful examination of these numerical values leads to the 
general rule. 


An 0 _ 
-^1 0" 


{377—5} 


(77=2,3, . . . °o) 


7i!4" 

and from the boundary conditions, equations (11), 

(77=2,3, ...«>) 


{377—5 } 
777714"”^ 


An 0 
-^0 0 


where by definition 

{377— 5} = 1X4X7X10X13X • • ■ X(3t7-5) 

In the expression for the local Mach number distribution 


6 


KEPORT 1149 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


evaluated at the crest of the wavy wall (a:=0, y=0), there 
occurs the following power series: 


F^'^nA2 g k" 

71 = 1 


This power series can be expressed in the closed form 


F=2 — 2 




The graph of F against t is a semicubical parabola with the 
cusp point at k=^ and F=2. With the necessary condition 

U 

that one and only one value of k correspond to a given value 
of F, the transonic similarity parameter k cannot be greater 
than 4/3. Moreover, the lower limit of k is zero when the 
amplitude of the wavy wall is zero but is different from 
unity. 

KECUHSION FORMULA FOR A 

With p = 1 , the recursion formula for A" / is (compare 
with eq. (69) of ref. 1) 

4(n+l)^V-2(7i+2)(g + l)5;^J+l+(g + l)(g + 2)5»-J-M"+J= 


5,=0 m=s-j, 


where 


(t?,= 1, 2, . . . 9=0,1, . . . n) (33) 


(g=n-l or n) 
“ 1 {q7^n—l or n) 


0 (q = n) 

1 iqT^n) 


From this relation, by equating coefficients of equal powers 
of k on both sides of the equation, the following equations 
are obtained: 

An-l 1_ £ 

8(n— 2)!4”-' 

In -1 1_ ^ ^ A\0 i_ ^ A1 0 

8(71-2)14"-* 8(77-3)14"-* ”* 


1 77" 

8 (77-3)14" 

1 77"-" 

16 (t7 — 4)14 


~,\Ail+3{Air)^]- 


(Al ?r 


(t 7=2, 3, . . . oo)_ 


Consider now g=77— 1; equation (33) becomes 

4(77+1).t4;;_5=277(77+ 2)A\ * — 

^77y)(777 + l)(77 — 777 — D^^t,* ° M"„-"„-_*i* — 
^ m=0 


^ n(n+l)AV A\y\° {n=l,2, . . . ~) 

Multiply both sides of this equation by ?■” and sum from 
77 = 1 to 77 = CO . Then 


4^ (?7-l- 1)^; _} 7-" = 2~^, 77 ( 77 + 2)A\ * r"- 

n=l 71 = 1 


;y^nr;^(m + l)(n-m—l )A’"+„' “ A ’‘„-"„-_\^ — 

■"n=l 77i=0 

■i ^*„“S l)A"+}i'> r" 

^ 7i = l 


The solution of this recursion formula proceeds as in the 
case p=0, the solution starting with 9=77 and descending 
towards g=0. Thus, for q=n, equation (33) becomes 


^V=-;|77^*„«^"+r (77=1, 2, . . .co) 


^V-i'=-g(«--l)-4V AU, (77=2, 3, . . .co) 

Hence, inserting the expression for given by equation 
(19) gives 


Aox = JZA\^r A 

71=1 71=1 


A — ’’s"' An ^ A A’>^ 0 -y.n 

-^1 0 Z-i -^n -1 • ^2 0 X-« -^n -2 • 

71 = 1 71 = 2 

Then it can easily be shown that the second term on the 
right-hand side of equation (36) is — ^ rir'^A^ /H, „')' and 
that equation (36) can be replaced by 

4(rH. d'= 27 -Mo i"-t-67-Ho/-|r(7-M„ rA\° A,Y' 


8 (77-2)14 


^(H*o“)"+* (77=2, 3, . . . co) (34) 


§ a; :l ,i. 7 -— i (S *■ 


i=-g H*o“Zl(7i-l)H"_?r”-* 

O 71 = 1 

or, with the aid of equations (25), 

7-Ho>=-i A\°{rA, o'r 


I nil I I I I II II 
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Hence, from the use of equations (25) and equation (27) for A 2 0 ' it follows, after routine calculations, that the integral 
of equation (37) is 

I o')*+(Ci+4) [i (r A 0 ')^- 2r‘A, „"] j 

A2 O 

Ci==64- " 

Now from equation (18); 


where 


(Avr 

n” 


Then 


(rA, g. („_„,)|4— (AVr r- 


(38) 


S ^V-V 




I ” 

\ '* 

»-3 ' 1 A 




n-4 r" + 


4^3(»-3)!4"-" 16^^ (n-4)!4”-“ 

Equating coefficients of equal powers of r on both sides of this equation gives 

(n=2,3, ... 00 ) (40) 

From this equation the expressions for all the coefficients of the type AlzlAi with n=2, 3, ... 00 and r=0, 1, 
... 00 can be obtained. Thus, for example, ’ 


■(M-3)!'^(n-4)!“^(n-5)! 


G 


a„-i 1 n"-°r 20 

»-2i 32 2)!' 

1 7i** ~ ^ C 4 14 

A=^^=-^^j[(5n+l)A;+32A?] -^^^-^-l- [(5n-l-l)A;+32A?] [(37n-b5) A ?+ 


(41) 


For q—n—2, equation (33) becomes 


4(n-f-l) A-2=2(a— l)(n-f-2) A-1 — n(n— l).^*,'— [(w + 2)(w-~m-2)^”+J* " ^V-’V-V + 

^ m=0 

(TO-|-l)(n-m-l)^“+J » A’‘-r'-l2^]—~n{n+l)A\° 24”+_V 


(42) 


The solution of this recursion formula follows along the same lines as that for equation (36) and leads to the following result: 

^n-. 1 (r 162 414 205 47 1 , .^sf 21 , 24 ,9 1 Hn"-*, 

512(L(n-3)!’^(n-4)!'^(n-5)!+2 + [(^^^+(^^=^ + 2 (^^J 


_(n-3)!”^(n-4)!'^(n-5)!^ 2 (n 


__J G "1 C 0N»+1 

,(n-4)!^(n-5)U 4"-^^ 

(71=3,4, ... 00 ) 


(43) 


From this equation there follow, in the usual manner, expressions for coefficients of the type A"Zlr+i- Thus, for 
example, with r=0, 

L\64 16 6 / (71— 3)!^V 64 "*"2 +48 / (n— 4)!+(,128‘*"32 ^ "^16 ^ j (71— 5)!"*~256 (71— 6)!j 4"-> 


( 44 ) 
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RECURSION FORMULA FOR 

With p=2, the recursion formula for is 
8(n+2M»/-2(n+4)(g+l)5'*tM;+?+ 

(?+!)(?+ 2)5" 

1 9 n+1-91 

= — 7 rn 3 *+=^ 5»-‘’,'‘’"+' 2 (m — 2 )( n— m + 

ffi >= 0 Tn = j — 

2)A" -“+“ “ i ri5" + V K "i(«— 

mM" ^ n(n+ 1)(5" A’‘ V ^+ 

gn+l,«+2 _4y ^n+l 0_^g0.n+2 

w(w + 2)(5%+= ay "+5“ AY ^%+_V) 

(n=l,2,. . .oo; g = 0,l,. . .71 + 2) (45) 

where 5 is defined in the usual manner. The solution of 
this recursion formula proceeds as in the previous cases, 
starting with g=7i+2 and descending towards 2=0. Thus, 
for g=7i,+2, equation (45) becomes: 

Aul=-^nAVA\%Y (71=1,2,. . .a.) 

or 

A'‘Y^=-~(n-2)AYA"A (71=3,4,. . .oo) 

Then, from equation (19), 

2) ^ (Avr - 

(71=3,4, ... <x.) (46) 

or 

An-2 2 !.2r 2)77 / .jg J.2r'Y^ 

By equating coefficients of equal powers of k on both sides 
of this equation, the following equations are obtained: 


A": 


(71 — 2)77" ■ 


An -2 2_ 
n 3" 


An -2 2 _ 
n 4" 


^!4n +2 

(n^— 4)n" 
n!4”+^ 

(n^— 4 )ti” 


/I 1 0 
^0 1 


7714" ■ 


(77+i)(+j;)*J 


^ (47) 


Consider now 3=77+!; equation (45) becomes 
8(77 + 2)+; +?= 2(77 + 2) (77 + 4)+; +i- 

i 77 S (»W+ 1)(77 — 777— l)+’»+;^'>+V-”+fl“- 

Z m—0 


1 


77(77+ 1)(+V+V+' + 24\‘+"+’ “)- 


77(77 + 2)(+V+“/+^“++^‘’+" Y ”) 

(77=1,2,. .. 00) (48) 


With the introduction of 


+ -22 2 j-^"+2 +2 0 2!j-^n-2 


n=2 


A >y ' An 2 A An 1 7 .* 

-^-1 2 .^L-C^n+l ' ^0 1 2: n ' 


+io=Z)-4;-??-" 


equation (48) is replaced by the following ordinary linear 
differential equation: 

8(r°+_i 2)'=2r“(r+_2 2')' + 12r(7’+_2 2')+16^^-2 2— 

lr\AA, „'+_2 20-^^'oV+o /'-i+\VM, 


+V 7-+: o " + ^ V^i o'-+V r +2 o " + 

+ 20 + ' + 10+20 

jri. I 212 0 ^ 0^0 


(49) 


By repeated use of equations (25) , the solution of this differ- 
ential equation is found to be as follows: 

8rM_, 2 =(+w{^ (r+, o')-~ (^+1 (r-+x 0 ')^- 


768 


2048 ” ^^”*”l6 “ "*”8 

^ (Ci + 4) [I (r+, „0^-++, o"]| i-(^^T+ 

+V [+. o-| (^+1 o')^]-»’+V+V (50) 

Finally, with the aid of equation (38) and the definition of 
+-I 2 , 


An -2 2 j f 

(256 


~ 1 1 2 "1 77"-" 

_(a-l)! (77-2)!'*‘(77-3)!j 4"-i‘ 


1 77" i_ , , r 1 

64 (a-4)!4"-*”'"512 L(«-l)! 

~| Vl \ O ')™ +2 

(w-2)! J 4"-+ 

("77 = 3, 4, ... 00; Ci = 64 


(^ 


,2 0 N 

From this equation, there follows in the usual manner 
formulas for the coefficients of the type +;rf+ 2- Thus, 
for r=0, 

A . -22 r 1 1 2 ~l 77"-^ 77 "-° 

"■‘'“L(ii-l)! (a-2)!’^(a-3)!j4"+^ (t7-4)!4"+^^ 

(^=3,4,... «.) 


DISCUSSION OF CONVERGENCE OF SMOOTH TYPE OF 
POTENTIAL FLOW PAST WAVY WALL 

In the preceding sections, a number of examples of recur- 
sion formulas and their solutions have been given in consider- 
able detail. The purpose of this exposition is threefold: 
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First, to show the inherent elegance of the method of integra- 
tion in series although the equations concerned are nonlinear 
in character; second, to present a type of analysis which may 
be useful in other problems involving nonlinear differential 
equations; and third, to indicate that an analytical proof of 
convergence may ultimately be obtained by careful examina- 
tion of the recursion formulas for the quantities and 

their solutions-- One example is the obtaining of the general 
expression for the coefficients A" “ and the subsequent con- 


clusion that k < - 


In actual practice it has been found more convenient to 
evaluate the coefficients A’^ f from their separate recursion 
formulas rather than to derive the general formulas. The 
appendix contains the exact numerical values of the coeffi- 
cients necessary for the development of the functions /„ to 
the eighth power in the transonic similarity parameter k. 
These coefficients are utilized to demonstrate numerically 
the test for convergence of the smooth symmetrical type of 
potential flow (eq. (3)) assumed in this report. Thus, when 
the form of f„ given by equation (8) is inserted into equation 
(5) for the local Mach number distribution and evaluated at 
the surface of the sinusoidal wall (y=0), the following result 
is obtained: 


M ^-1 

l—MJ 


» m 

l+Zl k” S (n- 

71=1 771 = 0 


-2m) cos (n—2m)x'^ A" ' 

p=Q 


(52) 


where 



denotes the integral part of n/2. 


At the crest of 


the wavy wall (a;=0), the point of maximum fluid velocity, 
with the numerical values for the coefficients inserted, equa- 
tion (52) becomes: 


M^-i 


' + S + + ‘‘+6W^256 ‘‘+ 


359381 ,, 7326757 ,, 

270X256" ‘*‘6480X256" ^ 


81688733 
”86400X256" ^ 

The critical value of k (that is, when M=\) calculated from 
this equation is 

A:„ = 0. 83244 


(Note that in ref. 1 the value A:cr=0. 83770 corresponds to the 
first six terms of eq. (53).) Consider now the infinite series 



where 


-/l-i 0 ^ 

71 = 1 


(54) 


An 0_ {_§^L__AI 


(n—2, 3, . . . co) 


An 0 
•^0 0 


The Cauchy ratio test o yields in the limit 71 — >oo 

•^0 0 

the result that the radius of convergence Ri is equal to or 
less than 4/3. If the corresponding ratios /? 2 n are formed for 

R 

the right-hand side of equation (53) and the quotient iV„ 
is calculated, the resulting sequence of numbers is as follows: 


n 


R2n 

Nn 


16.0 

8.0 

2.0 

2 

4.5 

1. 92 

2. 34375 

3 

3. 04762 

1. 78041 

1.71175 

4 

2.5 

1. 33366 

1. 87472 

5 

2.21546 

1. 35000 

1.61104 

6 

2.04167 

1. 17721 

1. 73132 

7 

1. 92481 

1. 19588 

1. 60950 


The noteworthy feature of this table is that, although 
Rin (and presumably /? 2 n) is converging quite slowly toward 
4 

/fi = 2 (and R 2 ), the quotient N„ exhibits a strong tendency 

to approach an asymptotic value N for a relatively low value 
of n. Figure 2 shows this tendency in a graphic manner. 
The apparent asymptote represented by the straight line is 
the ratio of 4/3, the limit of Ri„ as , and of 0.83244, the 
critical value of k. Certainly, the rapid approach of the 
lower dotted curve toward the apparent asymptote and the 
decreasing oscillations represented by the upper dotted curve 
indicate that the critical value of k is the radius of conver- 
gence Ri of the power series on the right-hand side of equation 
(53). Thus, the critical value of the stream Mach number 
marks the limit of convergence of the smooth symmetrical 
type of potential flow assumed. The ability to approximate 
closely the limiting value is a matter of luck; namely, the 
choice of the known comparison series. Once, however, 
the proper comparison series has been selected and the 
approach to an asymptote indicated, there can be no question 
of the meaningfulness of the approximate value of N obtained. 
It may be that one would like to extend figure 2 to n = 9. 
This extension would entail the forbidding calculation of an 
additional 185 coefficients A" ?. The result presumably 
would be to decrease slightly the approximate critical value 
of k and thereby raise slightly the straight-line asymptote of 
figure 2. This extension of figure 2 would show still more 
convincingly the approach to an asymptote and the dying-out 
of the oscillations. Perhaps more important still, figure 2 
definitely shows that conclusions based on less than six or 
eight terms are mere speculations in this field. 

CONCLUDING REMARKS 

If the numerical test of convergence presented is accept- 
able, the conclusion to be drawn is that smooth symmetrical 
potential flow past the wavy waU exists only for the purely 
subsonic range. Moreover no such flow can possibly repre- 
sent the transonic or mixed type for which a local region of 
supersonic flow near the solid boundary is imbedded in the 


■■I I I 
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otherwise subsonic stream. It follows as a corollary that 
the transonic or mixed type of flow past the wavy wall is 
necessarily an asymmetric one. This asymmetry in the 
flow pattern entails a resistance usually defined as wave 
drag. As shown by experimental observations, the shock 
wave associated with wave drag closes the downstream 
portion of the local supersonic zone. 

As a final remark — although the analysis and conclusions 
of the present work refer directly to the wavy wall, the 


suggested residt that the critical stream Mach number 
marks the limit of smooth potential flow very likely applies 
to other boundaries. This conclusion is based on the idea 
that the gradual transition from a purely sinusoidal wall to 
a boundary composed of a single bump, say, introduces no 
essential changes in the analysis. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 6, 1952. 


APPENDIX 


Coefficients for/u 


A 


1 

0 


o = il 

' 256 


A1 0 _ 
.^0 2 - 


1861 


' 3 X 256 ' 


A1 0_ 
-^0 3“ 


4896755 

'81X256® 


Coefficients for/j: 

A2 0 

^0 0 


A 


2 

0 


0 

1 


419 

72X256 


Al 


234215 

432X256' 


Al 


„ 42791533 

'■“648X256® 


Al 


0 

0 


1 

8 


Al 


0 

1 — 


11 

1024 


Al 


0 

2 


4085 

24X256' 


Al 


21287 

'■“324X256' 


NUMERICAL VALUES OF THE COEFFICIENTS A;’ 


All 

5 

Al 2_ 

65 

Al 1 = 

3385 

-^0 1 — 

256 

^0 2 “ 

'576X256 

■^0 3 

“1152X256- 

All 

2765 

Al 2_ 

26435 

Al 3_ 

3907 

^0 2 — 

9X256' 

3“ 

'81X256® 

^1 3 

'432X256' 

Al 1 

184345 

Al 2_ 

7 

Al 3_ 

23 

Jify 3 

1728X256' 

■^1 2 — 

72X256 

*^2 3 

'12X256' 

/4 I 1 

1 

Al 2_ 

53995 

Al 3_ 

583 

A\ 1 = 

64 

^1 3“ 

1728X256' 

-CI 3 3 — 

54X256= 

All 

33 

Al 2_ 

9 

Al 3_ 

119 

A-i 2 — 

64X256 

2 — 

32X256 

■^4 3 

12X256= 

All 

139 

Al 2_ 

2765 

Al 3_ 

1 

■^1 3 — 

4X256' 

■^2 3” 

96X256' 

•^5 3 

96X256 


8X256 


8X256' 


Al 


A2 1 
^1 3 


125 

“ 36X256 

A2 2 
-^0 2 

45895 

108X256' 

A2 2. 
-^0 3' 

4210409 

81X256® 

A2 2. 
-^1 2' 

5 

“ 256 

A2 2 
■^1 3 

815 

576X256 

A2 2 
2 

456185 

3456X256' 

A2 2, 
-^2 3 

1 

“ 128 

A2 2 
^3 2 

11 

32X256 

A2 2 
^3 3 

4811 

'“ 192X256' 

A2 2 
2 


A2 2 
-^4 3 


12245 

144X256' 

35023583 

1620X256® 

6245 

72X256' 

104063 

90X256® 

79 

3X256' 

167 

4X256' 

47 

192X256 

24199 

864X256' 

1 

12X256 

11 


A2 3 
.^0 3 


17999 

90X256® 


72577 

5760X256' 


45409 

2160X256' 


A 


2 3 

3 3 


3515 

864X256' 


A2 3 
-^4 3 


10229 

864X256' 


Ai\ 


73 

8X256' 


-^6 3" 


25 

12X256' 


2X256- 


? 
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Coeflicients for /a: 




A3 0 ^ 

A) 0 — 72 

., , 1765 

^ * 3X256" 

.3 a 475823 
“ " 15X256" 


.an 23603 
27X256" 

, 1351393 

2880X256" 

.3 3 _ 74353 
320X256" 


9843883 

81X256" 

.a . 155 

* 32X256 

.3 2 3083 

"" 24X256" 


AS 0 ^ 

^10-24 

, 30365 

" 48X256" 

J 3 2_. . 

2X256" 


.an 259 

72X256 

.3 4_ 117 

" ‘ 32X256 

.3 2^ 1785 

40X256" 


.an 149317 

432X256" 

29545 

96X256" 

. 3 2 _ 25 

2X256" 


. 3 0 1 

^20 32 

A3 1 1 

" " 256 



. a 0 __ 33 

32X256 

.an 139 

"" 2X256" 

A3 1 55 

" " 256" 



Coefficients for /. : 




7 

A 4 0__ * 

““ 6X256 

35251 

“ ‘ 90X256" 

.42 741115 

5184X256" 


390547 

720X256" 

80872738 

675X256" 

38381086 

405X256" 


1704729613 

16200X266" 

. 4 , 2465 

‘ ' 576X256 

423833 

1152X256" 


7 

.4 0 ' 

‘ 0 384 

30596538 

135X256" 

^V2 20845 

"^""■"144X256" 


32947 

' ^ 36X256" 

. 4 , 455 

" 96X256 

.42 _ 173 

12X256" 


20527213 
' " 162X256" 

, 127825 

" " 192X256" 

.4 2 30 

As 2 256" 


Ai 0 3 

"® 128 

.4 3_ 119 

" * 48X256 

A* 2 36 

®" 5X256" 


.4n 617 

" ‘ 288X256 

51343 

" 216X256" 



369587 

" " 1728X256" 

. 4 , 25 

" 48X256 



.4 0 _L 

^3 0 gg 

..n 11 

" " 24X256 

..n . 4811 

"" 144X266" 

.4,_ 275 

*" 8X256" 





12 

REPORT 1 149- 

Coefficients for/s: 


7 

A5 0 ‘ 

“ “ 3840 

AG 1. 
^0 1' 

90354659 
' 259200X256' 

Ar^ 

A ® ^ 

' 0 76g 

A5 1. 
■^2 1 

886243 

' ' 1080X256' 

A5 1 
■^3 1 

j 5 o_ 25 
' “ 1536 

A-21 

24375 

32X256' 

A ® ^ 

^3 0 334 

8950 

" 27X256' 

.3 0 25 

24X256 

.50 1375 

* ' 24X256' 

All 

Coefficients for /e : 


432X256 

AG 1 

•^0 1 — 1 

„ 285522407 

^ ' 4860X256' 

AU=- 

.0 0 91 

72X256 

1 

JL 

144862403 
' ' 810X256' 

AG 1 

-^3 1 

.0 0 23 

8X256 

A6 1 

-^4 1 

2553091 

" 2880X256' 

AG 1 

/I 5 1 

^oo_J^ 
' " 1024 

AG 1 

1 


n 56387 
" ' 96X256' 


AG 0 

^4 0 

AG 0 

-^4 1 — 


AG 0 _ 
■^6 0 


8X256 

815 

4X256' 

9 

20X256 
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I Coefficients for/?: 


5470394 

'81X256=' 

88735 

'96X256' 

25X15587 

'288X256'' 

25X997 

"24X256' 


'64X256 

__9 

"32X256 


2483629 

54X256' 

4371541 

5760X256' 

265901 

192X256' 

403535 

288X256' 

79841 

96X256' 

2205 

8X256' 

2401 

60X256' 


Coefficients for /g: 


A7 0 

■^0 0 — 

A7 0 

-^1 0 — 

A7 0 

^2 0 — 

A7 0 

-^3 0 


126X256 


'18X256 

287 

"144X256 

833 

"288X256 

343 

"144X256 


A7 0 

n — 


320X256 


A7 0 

-^6 0 


^45X256' 


13X19 

"l8X256' 

13X19 

"576X256 

89 

"64X256 

59 

“24X256 

47 

"18X256 
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